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ABSTRACT 

VLBI observations have shown that the parsec-jet of 3C 345 is formed by several 
components, ejected from the core with superluminal velocities and travelling along 
bent trajectories on the plane of the sky. We interpret the differences in velocity and 
position angle among the different features at formation time as the result of parsec-scale 
precession of the relativistic jet and calculate the aperture angle of the precession cone, 
the angle between the cone axis and the line of sight and the Lorentz factor associated 
with the jet bulk motion. We assumed a precession period of 10.1 yr, which is one of 
the -B-band light curve long-term periods reported in the literature. We propose that 
boosting of the underlying jet emission, which is time-dependent due to precession, is 
responsible for this long-term optical variability. Jet precession with periods of several 
years can be produced in super-massive black hole binary systems, when the secondary 
black hole is in an orbit non-coplanar with the primary accretion disk, inducing torques 
in the inner parts of the disk. Assuming that this mechanism is responsible for the 
jet precession in 3C 345, we estimate upper and lower limits for the masses of the two 
black holes, as well as their mean separation. We found a correlation between the 
formation of jet components and the occurrence of strong optical flares, as well as a 
very strong anti-correlation between the intensity of these flares and the time required 
for the components to reach the maximum flux density at radio frequencies. 

Subject headings: galaxies: active — galaxies: individual: (3C 345) — galaxies: jets — 
radio continuum: galaxies 



1. Introduction 

The quasar 3C345 (z=0.5928; Marziani et al. 1996), also known as 1641-^399 or 4C 39.48, was 
optically identified by Goldsmith &; Kinman (1965), and it is now known to be the nucleus of an 
elliptical E3 galaxy (Kirhakos et al. 1999). Strong variability (typical of an OVV object) is found 
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in its optical light curve, and periodicities of about 5 and 11 yr have been reported (Webb et al. 
1988; Kidger 1989; Zhang, Xie & Bai 1998; Zhang et al. 2000), although it has also been suggested 
that this variability may originate from non-linear and non-stationary stochastic processes (Vio et 
al. 1991). 

The radio flux density has been monitored at several frequencies (Waltman et al. 1991; AUcr, 
Aller &; Hughes 1996; Terasranta et al. 1998). The continuum spectrum is flat up to 10 GHz and gets 
steeper towards higher frequencies, with spectral index ranging from -0.9 (between 10^^-^ and 10^^-^ 
Hz) to -1.4 (between 10^^-^ and 10^^-^ Hz) (Bregman et al. 1986). Outbursts have been detected 
from centimeter to infrared wavelengths (Bregman et al. 1986; Stevens et al. 1996), besides time 
variable linear polarization at radio and optical frequencies (Moore & Stockman 1981; Bregman 
et al. 1986; Brown, Roberts &: Wardle 1994). In the optical regime, the degree of polarization 
decreases monotonically toward shorter wavelengths (Sitko, Schmidt & Stein 1985; Smith et al. 
1986; Mead et al. 1988; de Diedgo et al. 1994) and in the X-ray range, the source is weak and 
possibly variable (Halpern 1982; Makino 1989; Worrall & Wilkes 1990). 

VLA observations showed a faint halo around a bright core and an extended kiloparsec-jet 
(Kollgaard, Wardle & Roberts 1989). At parsec-scales, 3C345 exhibits a stationary core (Bartel et 
al. 1986; Tang et al. 1990) and a jet with superluminal components, which travel apparently along 
curved paths with variable velocities (Unwin et al. 1983; Biretta, Moore & Cohen 1986; Zensus, 
Cohen & Unwin 1995; Ros, Zensus & Lobanov 2000). These characteristics had been interpreted as 
a result of helical motion of the components along the jet, either in a pure phenomenological model 
(Qian et al. 1996) or as a consequence of HD or MHD instabilities (e.g., Konigl & Choudhuri 1985; 
Camenzind 1986; Hardee 1987; Qian et al. 1991, 1992; Camenzind & Krockenberger 1992; Steffen 
et al. 1995; Hardee 2000). 

An alternative approach is to assume that the helicoidal appearance of the jet is a consequence 
of precession. In this case, and considering that the bulk velocity is very large, each plasma element 
in the jet would move in an almost straight trajectory, defined by the jet direction when this element 
was ejected. From the observational point of view, superluminal components can be interpreted 
as these plasma elements, even if they are the results of shock waves propagating along the jet; 
their velocity should not be very different from the bulk velocity (e.g., Marscher & Gear 1985) and, 
at least close to their formation epoch, they should reflect the precessing jet direction. Curved 
trajectories, as found by Zensus, Cohen k, Unwin (1995) and Lobanov (1996) in 3C345, can be due 
to several reasons: jet instabilities (e.g., Hardee 2000), influence of the external environment (one 
of the mechanisms proposed by Walker, Benson & Unwin 1987 to explain the bending of the jet in 
3C 120), or even the superposition with other components formed at different epochs and moving 
with different superluminal velocities. In any case, the parameters of the precession model can be 
completely determined by the superluminal velocities and position angles of the jet components 
close to their formation epoch. 

The parsec-jet precession scenario has already been applied to 3C279, 3C273 and OJ287 
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( Abraham 8z Carrara 1998; Abraham & Romero 1999; Abraham 2000), where the velocity variations 

among components were attributed to differences in the angle between their trajectories and the 
line of sight. Besides, it was found that the periodic outbursts seen in the optical light curves of 
3C 279 and O J 287 could be due to changes in the boosting parameter (Abraham & Carrara 1998; 
Abraham 2000). Although this model explained well the behavior of these objects, the precession 
parameters were obtained from a limited set of data and should probably be revised to include 
new and more precise data (Homan et al. 2001; Wchric et al. 2001). We will show in this paper 
that the variations in the ejection angles and velocities of the jet components of 3C 345 can also 
be interpreted as the result of parscc-jct precession. Furthermore, we show that the presence of 
long-term optical variability can be attributed to variable boosting of the non-thermal radiation 
from the underlying jet. 

In Sect. 2, we describe the precession model used in this work. In Sect. 3, we present the 
precession model parameters for 3C 345, as well as the influence of precession in the absolute core 
position shifts due to opacity effects. In Sect.4, we discuss, in the framework of the precession 
model, other observed characteristics, such as the relation between optical flares and the origin 
of supcrluminal components, 5-band long-term periodic variability, polarimetric observations and 
dispersion in the apparent angular size of the components. Furthermore, orbital parameters of a 
possible super-massive black hole binary system in the nucleus of 3C 345 are estimated based on 
the precession model, with additional constrains from optical continuum and spectroscopic data. 
Finally, conclusions are presented in Sect. 5. 



Let us consider a relativistic jet with bulk velocity /3, precessing with a constant angular 
velocity cu and period P around an axis, forming a conical surface with aperture angle Q. The cone 
axis forms an angle (pQ with the line of sight and presents a projected angle r]Q on the plane of the 
sky. The instantaneous position of the jet is represented by the angles (p and r]. Their dependence 
upon the time t', measured in the comoving frame, can be expressed as (Abraham &: Romero 1999; 
Abraham 2000): 



2. 



Precession model 



r]{t') = arctan(— ) 



(1) 



= arcsmi 



(Vx^Ty^) 



(2) 



with 



X = (cos O sin (po + sin CI cos sin ut') cos rjo — sin J7 cos ut' sin rjo 



(3) 
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y = (cos sin 0o + sin cos <^o sin wt') sin r/o + sin O cos oot' cos % (4) 

The apparent velocity /3obs (i^i units of light speed c) is related to the proper motion of the jet 
components through: 



where z is the redshift and Dj^ is the luminosity distance, defined as (Carroll, Press &: Turner 1992): 



Dl = ^^^^^(f^M,OA,z) (6) 
-no 

where Hq is the Hubble constant and £^ is a function that depends on z and the dimensionless 
density parameters Om and Oa- Throughout the paper, it will be adopted Hq = 71 km s~^ Mpc~^, 
= 0.27 and = 0.73, as derived from the recent PFMAP results (Bennett et al. 2003). In this 
case, since J^m + f^A = 1 the function E is given by (Carroll, Press & Turner 1992): 



For the chosen cosmology, = 2A7h~^ Gpc, with h = i^o/lOO, from which 1 mas yr"^ 
corresponds to an apparent velocity of 2A.5h~^c. 

The viewing angle (f) is not measured directly, but it is related to (3ohs by- 



l-/3cos[(/.(t')] ^ ^ 

The elapsed time between two events in the observer and comoving frameworks {At and At' 
respectively) are related by the Doppler factor S as: 



with 
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^(</',7)=7"Hl-/3cos[0(O]} 



-1 



(10) 



and 



7 = (1 - 



(11) 



Thus, given J7, 0o, rjo, (3 and P, we could predict at any time the observed properties of a 
precessing jet. 

The flux density in the optically thin regime will be boosted according to the expression: 



where a is the spectral index {S^ (x v °), p = 2 for a continuous jet and p = 3 for discrete features 
(Blandford & Konigl 1979; Lind & Blandford 1985). 



In this work we used the VLBI data of 3C 345 (core-component distance r, position angle on 
the plane of the sky rj and flux density Si,) found in the literature, between 5 and 22 GHz, covering 
almost 20 yr of monitoring (Biretta, Moore &; Cohen 1986; Unwin & Wehrle 1992; Unwin et al. 1994; 
Brown, Roberts &; Wardle 1994; Zensus, Cohen & Unwin 1995; Leppanen, Zensus &; Diamond 1995; 
Lobanov 1996; Ros, Zensus & Lobanov 2000). We adopted the emergence epochs for the identified 
components given in the references, except for C7, for which we ruled out those observations with 
not clear identification (e.g., C76 in Unwin & Wchrlc 1992 and Unwin ct al. 1994). As there are 
no observations of the older components C2 and C3 close to their ejection epoch and C9 has only 
one reported observation (Ros, Zensus & Lobanov 2000), we did not use them in our fitting. In 
order to estimate the apparent proper motion /j, and position angles of the components, it was 
assumed quasi-ballistic motions in the inner region of the quasar (r < 1.0 mas), since the data 
show that the trajectories are bent beyond this distance. It is important to emphasize that some 
components, such as C4 and C5, present variable proper motions and/or position angles even for 
r < 1.0 mas (e.g., Lobanov 1996); in those cases, we have considered only observations taken before 
the occurrence of significant changes in those quantities. We should note that all the superluminal 
features became optically thin at distances smaller than 1 mas, since this fact is important for the 
discussion in Sect. 4.1. 



5j(z.) = 5j(^.)(5(</.,7)^+' 



a 



(12) 



3. Determination of the precession parameters 
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As we used data obtained at different frequencies, it was necessary to take into account opac- 
ity effects on the determination of the absohitc core position, which introduce frequency-dependent 
shifts in the core-component distances and proper motions. In the case of a processing jet, these 
corrections are time dependent and its magnitude calculated in Appendix A. The parameters in- 
volved in this calculation are the integrated synchrotron luminosity L^y^-, the ratio between upper 
and lower limits in the relativistic particle energy distribution 7max/7min, the intrinsic jet aperture 
angle V'', and a constant parameter k^- Following Lobanov (1998), we assumed 7max/7min = 100, 
fee = 1 and a constant value of 1.58 x 10^^ erg s~^ for Lsyn, even though this quantity could be time 
variable. Finally, we chose ifj' = 1° based on the observed angular size of the jet components as a 
function of their distances from the core. 

Among the parameters discussed in Sect. 2, we constrained the precession period P using the 
5-band photometric data. In fact, Zhang, Xie k. Bai (1998) reported a period of 10.1 ± 0.8 yr 
in the light curve of 3C345 after application of Jurkevich test (Jurkevich 1971). We assumed 
that this periodicity is due to variable boosting of the jet emission as the Doppler factor varies with 
precession (equation [12]). 

We adopted an unique Lorentz factor for the bulk motion of all components, compatible with 
the velocity of the fastest component 8.7h~^ c) given by jmin = (1 + /^obs)''^''^ (t ^ 12.3 for 
h = 0.71). After fixing 7 close to its lower limit, we selected the parameters CI, (f)Q and ijq that 
fitted the apparent velocities and position angles of the jet components using equations (l)-(4), 
with t' as implicit variable. Then, we checked the behaviour of hPohs and rj as functions of time 
through equations (l)-(5), assuming that the position angle of the jet component represents the 
jet position at the epoch when the component was formed. This procedure was repeated until a 
good fitting was obtained. For each set of precession parameters, we applied the core opacity shifts 
to the observational data, which introduced refinements in the kinematic and model parameters. 
The reduced chi-square value for both l3{t) and ij^t), for six supcrluminal components and three 
parameters is about 3, most of it due to uncertainties in the position angle ij. 

In Fig. 1, we present the difference in core position Arcorc between 5 and 22 GHz (the frequency 
range used in our work) as a function of time, calculated with the precession parameters. The mean 
value obtained for Arcore is 0.26 mas, while its lower and upper limits are respectively 0.21 and 0.31 
mas. Lobanov (1996) found an average offset of 0.328 ± 0.020 mas between 1992.45 and 1993.72, 
in good agreement with our estimate of 0.30 mas for the same interval. 

The kinematic parameters of the jet components corrected by opacity effects are presented in 
Table 1. The quoted errors correspond to the range for which a reasonable fitting for the data 
could still be found. 

The precession parameters that best fitted the data are given in Table 2, while Fig. 2 presents 

the model fitting in the (r/, h(3ohs)-> (t, h(3ahs), (t, f])-, {t, 5) and (t, ip) planes. Lower limits for 5 
(full triangles) calculated from X-ray observations (Unwin et al. 1983, 1997) are presented in the 
same figure. The Doppler factors predicted by our model are always above the lower limits imposed 
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by the X-ray observations. In fact, although other combinations of the precession parameters (f)o, 
r/o, ^ and 7 also fitted the kinematic data, they resulted in Doppler factors incompatible with the 
lower limits imposed by the observed X-ray fluxes. 

Another consequence of the precession model is a possible dispersion in the component angular 
sizes. As mentioned before, the observed sizes ip depend on the angle between the jet and the line 
of sight (equation [A2]). Observations show a dispersion 11° < V ^ 45° with a mean value of < 
il> >= 28°, at least for distances r < 1 mas, where we assumed ballistic motion. These observations, 
together with the modeled time behaviour of the apparent jet aperture angle (assuming an intrinsic 
aperture ip' = 1°) are presented in Fig. 2e. The precession model predicts 14° < V' < 42° and 
< ip >= 28° and, except for C3a, there is good agreement between the individual observations 
and the model, which seems to indicate that the jet components have similar intrinsic sizes. If 
we assume that the time behaviour of ip is fully described by precession effects, there are several 
alternatives to explain the discrepancy between the calculated and predicted values for C3a: its 
intrinsic size is larger than 1°, the estimated value of its apparent size is uncertain due to its low 
flux density or to its possible relation with C4 (Lobanov 1996), or finally because few observations 
are available for r < 1 mas. 

4. Discussion 

General results obtained from previous sections are discussed below. 

4.1. The optical light curve 

In Fig. 3, we show the -B-band light curve of 3C 345 at the observer's reference frame (Zhang, 
Xie &: Bai 1998). We also show, as a continuous line, the boosted emission of the underlying jet, 
calculated from equation (12) with aopt = 1-66 (Hagen-Thorn et al. 1996), p = 2, the Doppler factor 
obtained from the precession model and the flux density in the comoving reference frame Sj = 7 
nJy, which gives the right values for the observed flux density at epochs 1971-1972 and 1990-1991. 
The inferred value for the underlying jet flux density is an upper limit, since the total luminosity 
in the optical band is probably due to the superposition of other processes, such as emission from 
the accretion disk and from the superluminal components. The existence of a weak underlying jet 
was already postulated in 3C 279 and OJ287 (Abraham & Carrara 1998; Abraham 2000). 

From the light curve, and whenever observations were available, we were able to associate the 
emergence of jet components to the occurrence of flares in the optical band^, even though for the 
older jet components uncertainties in to lead to a non-unique association. The emergence time of 



^We have considered as optical flares short-time variations (between 5 and 90 days) that introduce changes in the 
flux density higher than 0.5 mjy. 
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the superluminal features and their associated flares are explicitly marked in Fig. 3. There we 
can also see that only flares stronger than a certain limit (about 1.5 mJy) produced detectable 
superluminal components. 

It is not clear whether these optical flares are related to short-lived stages of shock evolution 
(e.g. Marschcr & Gear 1985), or a consequence of some instability produced in the accretion disk 
(Romero et al. 2000). An unexpected relation between the flare flux density at optical wavelengths 
and the elapsed time between the epochs of superluminal component ejection and occurrence of 
maximum flux density at radio frequencies gave us some insight on the origin of these flares. 
Considering only the interval for which the movement of the jet components is quasi-ballistic, and 
deflning ATjnaxi^) as the the elapsed time between the epoch when the components reach their 
maximum intensity tj^axi^) at the observed radio frequency v and their formation epoch (to), we 
obtain: 



The values of AT^ax for C3, C3a, C4, C5 and C7, obtained from the flux densities given in 
Lobanov (1996), are presented in Table 3. Wc ignored C2 because the maximum occurred before 
the first observation, and C8 and C9 because few observations are available. We calculated ATmax 
at 22 GHz for all components except for C3 and C3a, for which we used data at 10.7 GHz, because 
they did not have good time coverage at 22 GHz. The error introduced by the use of a different 
frequency did not change appreciably our results. 

In Fig. 4a, we show the 5-band flux densities Sb of the flares associated with the jet components 
as a function of ATmax (i^)- This plot shows a clear and unexpected anti-correlation between these 
quantities (R = —0.97, where R is the correlation coefficient obtained from a power-law fitting). 
This kind of correlation was not seen when the maximum radio fiux density was considered. 

To verify whether this anti-correlation is real or a consequence of variable boosting, it is 
necessary to make the analysis in the comoving reference frame. As jet components are moving 
with relativistic velocities, the frequency of the emitted radiation in the comoving frame v' and the 
frequency u measured in the observer's reference frame are related hy u = 5v' . Besides, because time 
dilation, Arinax(i^) in the observer's framework is related to /STl^^{v'), defined in the comoving 
frame as: 



But what we want is a relation between AT^^x and ATmax always at same comoving frequency, 
for example at frequency v. Assuming that each component propagates with constant velocity in 
the comoving reference frame, Abraham (2001) showed that: 



ATmax (z^) = imax(i^) - ^0 



(13) 



^T^^{^') = AT^^{v/5) = 5ATmax(j^) 



(14) 
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Ar;,,(z./<5) = <5Vf AT4^(i.) (15) 

where the parameter b depends on the dominating energy loss process (Marscher & Gear 1985). 
For instance, for jet components in the adiabatic-loss stage, b is given by: 

^_ 2(2^ + 1) + 3a(8 + 2) 

3(s + 4) ^ 

s = 2ar + 1 (17) 

where ar is the spectral index at radio wavelengths, and the parameter a depends on the variation 
of the magnetic field with the distance y' in the jet's reference frame [B{y') oc y'^""]. For a = 1, 
B± « {B± and are respectively the component of the magnetic field oriented perpendicularly 
and parallel to the jet axis), while for a = 2, B\\ » B±. 

Substituting equation (15) in equation (14), we have that ATmax measured in the observer's 
reference frame and in the comoving frame at same frequency v are related by: 

Ar^a.(^) = d^'-'AT^^^ii.) (18) 

In Cols. 4 and 5 of Tabic 3, wc present respectively the Dopplcr factor So at the formation 
epoch to and the measured spectral index a^- Assuming that components are in the adiabatic phase 
and a = 1, we obtained b through equation (16) (Col. 6 of Table 3). With these parameters and 
equation (18), we determined ATl^^{i/), displayed in the Col. 7 of Table 3. 

The transformation of the optical flux density to its value in the the comoving reference frame 
depends on the nature of the flare. If it is associated with some kind of instability in the accretion 

disk, such as perturbations due to the passage of the secondary black hole through the disk (Romero 
et al. 2000), then the fiux density of the fiares would not be boosted. As it can be seen in Fig. 4b, 
the anti-correlation between optical fiux density and AT^^^^{iy) is maintained, indicating that the 
stronger the flare is, the faster the propagation of the perturbations in the disk is. 

Other possibility is that the optical flares are produced in the jet in the early stages of the 
shock evolution. In this case, relativistic corrections are needed and: 



5b = 5^+'^°^S^ 



(19) 
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where is the 5-band flux density of the flare measured in the comoving frame. The values of 
are shown in the last column of Table 3. 

In Fig. 4c, wc plotted S'-q in terms of AT^^^^v). Again, we adopted aopt = 1-66 (Hagen- 
Thorn ct al. 1996). Clearly the remarkable anti-correlation seen in Figs. 4a and b disappears 
completely, indicating that if the optical flares are short-lived phases of the evolution of superluminal 
components, the anti-correlation is only a consequence of boosting. However, if the optical flares 
originate from jet components, we would expect the anti-correlation between maximum flux density 
in the observer's reference frame and ATiaaxi^^) be observed also at radio frequencies. Since no 
correlation was found between ATmg^{i') and maximum flux density at 10.7 or 22 GHz, the scenario 
where the optical flares are produced in the accretion disk is favoured. It is important to emphasize 
that the results obtained in this section are not altered if we had considered a = 2 in our calculations. 



4.2. Polarimetric observations and precession model 

Polarimetric interferometry of 3C 345 (Brown, Roberts &i Wardle 1994; Leppanen, Zensus 
Sz Diamond 1995; Taylor 1998; Ros, Zensus &; Lobanov 2000) has shown that the electric field 
orientation and the fractional polarization vary along its jet. Brown, Roberts &; Wardle (1994), 
observing at 5 GHz, found that the magnetic field is almost aligned with the local jet direction, 
except at the position of the superluminal features. From observations at 22 GHz, Lobanov & 
Zensus (1996) and Ros, Zensus &: Lobanov (2000) proposed the existence of a transition region 
where the magnetic field changes its orientation from transverse to longitudinal. In the inner parts, 
strong shocks would dominate, favouring magnetic field perpendicular to the jet axis, while shocks 
in the outer region would be weaker, and they would not have too much infiuence on the magnetic 
field orientation. The increase of the fractional polarization with distance along the jet would be the 
result of the superposition, close to the core, of components with different electric field orientations 
(Ros, Zensus & Lobanov 2000). 

If there were alignment between the magnetic field and the underlying jet, as in the case of 
3C 345, jet inlet precession would introduce fast time variations in the polarization angle, specially 
when the viewing angle is small. In fact, at the epochs for which the viewing angle is smallest, our 
model predicts variations in position angle of 60 degrees during an interval as short as 2.6 yr (Fig. 
2c). The vectorial superposition of the electric fields associated with jet components formed during 
these epochs may lead to the cancellation of the total polarized flux and consequent decrease in 
the fractional polarization near the core. For larger distances, the differences in apparent velocities 
between the different components would decrease the superposition and the fractional polarization 
would be higher. 
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4.3. Physical parameters of a possible super-massive black hole binary system 

The precession period of 10.1 yr assumed in our model is certainly inconsistent with the period 
of 1000 yr calculated by Lu (1992) using the Lense-Thirring effect (Lense & Thirring 1918), for 
which jet precession is due to the misalignment between the angular momenta of the accretion disk 
and of a Kerr black hole. The precession of the accretion disk tidally induced by a secondary black 
hole in a black hole binary system seems to be a probable precession mechanism (Katz 1980, 1997; 
Romero ct al. 2000). Assuming that last interpretation is correct for 3C345, we can estimate the 
physical parameters of the binary system. 

Let us consider that the primary and secondary black holes, with masses Mp and Ms respec- 
tively, are separated by a distance rpg. From Kepler's third law, we can relate rps to the orbital 
period of the secondary around the primary black hole Pps'. 



3 _ GMtot 2 .20) 

where G is the gravitational constant and Mtot is the sum of the masses of the two black holes. 

In the observer's reference frame, the orbital period Pp^^ will be: 



Ppr = (1 + ^)Pps (21) 

If the orbit of the secondary is non-coplanar with the accretion disk, torques could be induced 
in the inner parts of the accretion disk, producing its precession. Considering that the outer radius 
of the prcccssing part of the disk is rj, Papaloizou &: Terquem (1995) and Larwood (1997) calculated 
its precession period in terms of the masses of the black holes: 



27r,, , 3/7-2n\GMs 

— (1 + 2;) = — ^=-^^^^cos6' (22) 

where n is the politropic index of the gas (e.g., n = 3/2 for a non-relativistic gas and n = 3 for the 
relativistic case) and is the inclination of the orbit of the secondary with respect to the plane of 
the disk. 

As it is believed that accretion disks are the main candidates for feeding jets, they are expected 
to form coupled systems (e.g., Donea & Biermann 1996, 2002). In this case, the jet would precess 
at the same rate than the disk (P = Pj), forming a precession cone with half-opening angle equal 
to the angle of orbit inclination (0 = 0). Thus, equation (22) becomes: 
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/7-2n^ 


G(Mtot - Mp) 











where Mg was replaced by Mtot — -Mp. 

Prom equation (23), can be obtained in terms of Mp and Mtot by: 



Stt / 5 - n \ (1 + z) r^^ 
'^'^ 3 V7-2ny Pcosi7 VGMt 

where Xp = Mp/Mtot- However, equation (24) is vaUd only if the disk precesses like a rigid-body, 
implying that must be appreciably smaller than Vps (Papaloizou & Terquem 1995). Thus, if Mtot 
and Pps are previously known, this condition can be used to put an upper limit to the mass of the 
primary black hole, as well as a lower limit to the secondary mass. 

Using the velocity dispersion obtained from the H/3 line-width and the optical continuum 
luminosity, Gu, Cao &; Jiang (2001) found a mass of 8 X 10^ Mq for the central object in 3C 345, 
which it is assumed to be the value of Mtot in this work. In Table 4, we present the separation 
between the two components of the binary system and upper (lower) limits for the mass of the 
primary (secondary) black hole for several orbital periods. We show in Fig. 5 the outer radius of 
the precession disk, as given by equation (24), as a function of Xp. 

We can rule out orbital periods P^^^ < 4 yr because, in that case, the time-scale for losses due 
to gravitational radiation is smaller than about 500 yr (e.g., Begelman, Blandford &; Rees 1980; 
Shapiro &; Teukolsky 1983), which would lead to significant changes in the orbit of the secondary 
and consequently to non-periodic disk precession. In addition, from Table 4 we should not consider 
periods P"^^ > 6.1 yr since they correspond to a mass for the secondary black hole larger than 
that of the primary, which seems to be physically unlikely. For these reasons, the possible ranges of 
primary and secondary masses are 4 x 10^ Mq < Mp < 5 x 10^ Mq and 3 x 10^ Mq < Mg < 4 x 10^ 
M0, respectively. 

According to Romero et al. (2000) , the passage of the secondary black hole through the primary 
accretion disk can produce density waves and induce shocks in the jet, leading to the formation of 
the superluminal components. Since the secondary black hole passes through the disk twice per 
orbit, the maximum efficiency is obtained when each passage produces at least one jet component. 
Assuming that only one component per passage is formed, the interval between successive ejec- 
tions of superluminal components, measured in the observer's reference frame P^^^, corresponds to 
P°^^/2. If not all the interactions are able to generate a jet component, the observed periodicity 
would be a multiple of half the orbital period {P°^^ = vPpg^/2, with v = 1,2,3, ...). From the 
emergence epochs of the jet components, we observe that the median ejection interval is 2.6 yr. 



,1/3 



(24) 
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Assuming that this interval corresponds to the maximum ejection efficiency, P^^^ = 5.2 yr, similar 
to the period of 5.1 yr found in the 14.5-GHz and 5-band light curves (Zhang, Xie & Bai 1998; 
Kelly et al. 2003). Moreover, if the 2.6-yr ejection period corresponds to twice the half period, 
we would obtain a full period of less than 4 yr, incompatible with the condition found for stable 
precession mentioned above. 

On the other hand, if each passage of the secondary object through the accretion disk generates 
multiple components Nc, the relation between P^^^ and P°g^ is P°^^^ = TP°^^/2, where T = v/Nc- 
Even in the simplest case where Nc = 2, and assuming again P^-^^ = 2.6 yr and f = 1, we obtained 
Pp^^ = 10.4 yr. However, this value for the orbital period corresponds again to a binary system 
where the secondary black hole is more massive than the primary, so we can also rule out this 
possibility. Therefore, it is reasonable to conclude that on average, each passage of the secondary 
through the accretion disk produces one superluminal component. 

5. Conclusions 

In this work, we used data available in the literature involving the parsec-scale radio structure, 
which comprise almost twenty years of monitoring of 3C 345. We assumed that all superluminal 
features have quasi-ballistic motions in the inner regions (r < 1 mas). 

We showed that, close to the core, the trajectories of the superluminal components can be 
explained quite well if the jet is precessing with a period of 10.1 yr, one of the B-band periodicities 
reported by Zhang, Xie & Bai (1998). As a consequence of this, we found a correlation between time 

variation of the Doppler factor and the optical light curve, indicating that probably the boosted 
optical emission from the underlying jet is responsible for the 5-band long-term periodic variability. 
An upper limit of 7 nJy for the intrinsic flux density of the underlying jet has been estimated in this 
work. Lower limits for the Doppler factor calculated previously from X-ray observations (Unwin et 
al. 1983, 1997) were used as constrains in our model. 

In this work, we also analyzed the influence of jet precession in the observed position of the 
optically thick core at different frequencies (Blandford Sz Konigl 1979; Lobanov 1998). We found 
that precession introduces a time modulation in the shift corrections, as well as in the direction in 
which the corrections should be applied. However, core position shifts calculated in this work did 
not exceed 0.31 mas. 

Considering that jet precession in 3C 345 is driven by a secondary super-massive black hole in 
a non-coplanar orbit around the primary accretion disk and using the total mass of the two black 
holes derived from the H/3 line- width and optical continuum luminosity (Gu, Cao & Jiang 2001), 
we estimated the masses of the primary and secondary black holes as 4 x 10^ Mq < Mp < 5 x 10^ 
M0 and 3 x 10^ M0 < Mg < 4 x 10^ Mq, respectively. In addition, we found that the distance 
Tps between them should be 5.5 x 10^^ cm < rpg < 7.3 x 10^^ cm and their orbital period 2.5 yr < 
Pps < 3.8 yr. 
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The observed dispersion in the apparent sizes of jet components could be attributed to jet 

precession and interactions between jet and environment matter. It is corroborated by sudden 
variations in the radio flux density of the knots when there are changes in their propagation direc- 
tion. Furthermore, those interactions seem to be responsible for gentle bending seen at distances 
larger than 1 mas, which could explain also the difference between jet orientation at parsec- and 
kilopar sec-scales . 

We were able to associate the emergence of jet components to the occurrence of flares in the 
optical band stronger than about 1.5 mJy. We also found an unexpected relation between the 
flare flux density at optical wavelengths and the elapsed time between the epochs of superluminal 
component ejection and occurrence of maximum flux density at radio frequencies, which led us to 
the conclusion that these optical flares are originated in the accretion disk and not in the shock 
waves seen as superluminal features. 

Finally, the complex polarimctric structure of 3C 345 in parsec-scalcs (Brown, Roberts & War- 
die 1994; Leppanen, Zensus & Diamond 1995; Taylor 1998; Ros, Zensus & Lobanov 2000) may 
be explained by our precession model, when the superposition of jet components with different 
position angles is considered. 

This work was supported by the Brazihan Agencies FAPESP (Proc. 99/10343-3), CNPq and 
FINEP. We would like to thank the anonymous referee for careful reading of the manuscript and 
for useful comments and suggestions. 



A. Opacity effects on core-component distance in a precessing jet 

The kinematic properties of the superluminal components in the jet of 3C 345 depend on 
the accurate determination of the core-component separation. As pointed out previously (e.g., 
Blandford & Konigl 1979; Lobanov 1996, 1998), the absolute core position rcore depends inversely 
on the frequency when the core is optically thick, which introduces a shift in the core-component 
separation. Following Blandford & Konigl (1979), we can write the absolute core position as: 



_ 4.56 X 10-^^(1 -Hz) 

^core(,^j — 7 ij^ 

Di^^^ke ip sin (j) 



L 



syn 



sin( 



J3{1 - (3 COS 0) ln(7max/7min) 



2/3 



V ^(mas) 



(Al) 



where z is the redshift, Z?l is the luminosity distance (in units of parsec), is a constant (fee < 1; 
Blandford & Konigl 1979), L^jn is the integrated synchrotron luminosity (in units of erg s~^), while 
7max and 7min are related respectively to the upper and lower limits of the energy distribution of 
the relativistic jet particles. The quantities -0 and v are respectively the observed aperture angle 
of the jet (in radians) and the frequency (in Hz); the former is related to the intrinsic jet aperture 
angle -0' through (e.g., Mutel et al. 1990): 
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tan('0/2) = tan('072) cot (j) 



(A2) 



The core position shift Arcore between frequencies ui and 1^2 (^2 > z^i) is given by: 



Arcore(z^l,I^2) 



4.56 X 10-12(1 + ^-j r gin ^ 



(mas) (A3) 



Di^j^kl^^i/j sin ^ L^(l - /3cos(^)ln(7i„a^/7mi„). 



Note that if wc substitute in equation (A3) a simpler version of equation (A2), ip ~ ip' csc9, 
we obtain equation (11) given in Lobanov (1998). 

We can see that equation (A3) depends on the angle between the jet and line of sight; in the 
case of a jet which is precessing, this angle is a function of time, what obviously introduces a time 
dependency in Arcore- On the other hand, the shifts in the core-component separations do not 
occur in a fixed direction, but they are oriented according to the direction at which the jet inlet 
is pointed. As the jet inlet is not resolved by observations, changes in its direction will reflect on 
changes in the position angle of the core region. Thus, a jet component, located at a distance r 
from the core and with a position angle rj, will have right ascension and declination offsets (Aa 
and A6 respectively) given by: 



where r/c is the position angle of the core at the epoch tobs in which observation is acquired. Once the 
precession model parameters are known, we are able to determine the second term of the equations 
(A4a) and (A4b) and correct the component position by core opacity effects; using them, we can 
determine the corrected core-component distance rcorr through: 



Note that if Arcore = 0, rcorr = r. Other particular case is found when there is alignment 
between the position angles of the core and of the jet component {rj = rjc), such that rcorr = 



Aa{ui, 1/2) = r{vi) sin[77(tobs)] - Arcore(z^i, 2^2) sin[r/c(tobs)] 
A(5(i/i, 1^2) = r{i'i) cos[r?(tobs)] - Arcore(z^i, 1^2) cos[77c(tobs)] 



(A4a) 
(A4b) 



{1^1) = V riui)"^ + Arcore(z^i, ^'2)2 - 2r(z/i)Arcore(i'i, 2^2) cos(?? - r?c) 



(A5) 



r — Ar, 



core- 
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Fig. 1. — Time behaviour due to jet precession 
of the difference between core position obtained 
at 5 and 22 GHz. 
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Fig. 2. — Precession model applied to the parsec-scale jet of 3C 345 with model parameters hsted in 
Table 2. (a) - (e): Continuous lines are model predictions on the planes (rj, h/3ohs), (t, /i/3obs)) (t, 
rj), {t, 6) and {t, ijj). Full circles and triangles represent, respectively, observations and lower limits 
for the Doppler boosting factors, calculated from X-ray observations (Unwin et al. 1983, 1997). 
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Fig. 3.— 5-band light curve of 3C 345 (full cir- 
cles). Solid line represents the expected contri- 
bution of the underlying jet. The uncertainty 
in the emergence epoch of each jet component 
is represented by the horizontal lines under the 
component label (see Table 1), while the verti- 
cal lines mark the flares that could be associated 
with these components. 
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Fig. 4. — Comparison between radio and optical properties of the superluminal components and 
flares. Filled symbols represent jet components for which 22 GHz radio observations were used, 
while open symbols correspond to 10.7 GHz observations, (a): 5-band flux density of the flares 
associated with the superluminal features as a function of the elapsed time between the epochs of 
component ejection and occurrence of maximum flux density at radio frequencies, (b): The same 
as (a), but as a function of the elapsed time in the comoving jet reference frame, (c): The same 
as (b), but considering boosting of the optical flux density. The solid line in each plot corresponds 
to the minimum least squares fitting using power laws. 
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Fig. 5. — Outer radius of the precessing disk 
function of Mp/Mtot- Each fuh hne corresponds 
to a different orbital period, increasing from right 
to left (see Table 4). The big circles show the 
values at which the outer radius of the precessing 
disk equals the separation between the two black 
holes. The dashed line marks the position in 
which the masses of the primary and secondary 
black holes are equal. 



Table 1. Parameters of the superluminal components of 3C345. 



Component 


to (yr) 


/i (mas yr ^) 




V n 


OF 


C3a 


1977.5 ± 1.5 


0.355 ± 0.025 


8.7 ± 0.4 


-98 ± 6 


yes 


C4 


1980.1 ± 1.3 


0.225 ± 0.015 


5.5 ± 0.3 


-120 ± 8 


ycs^ 


C5 


1982.3 ± 1.0 


0.200 ± 0.030 


4.9 it 0.6 


-107 ± 8 


yes 


C6 


1987.5 ± 1.0 


0.347 ± 0.030 


8.5 ± 0.6 


-98 ± 8 


no*" 


C7 


1991.0 ± 1.0 


0.208 ± 0.025 


5.1 ± 0.5 


-97 ± 6 


yes 


C8 


1993.0 ± 1.0 


0.196 ± 0.030 


4.8 ± 0.6 


-110 ± 6 


yes* 



Note. — to is the formation epoch, fj, and hfioha a^e, respectively, the proper 
motion and the apparent velocity in units of light speed c, and ?? is the ejection 
angle. Col. OF indicates whether or not there is some optical flare observed in 
the i?-band associated to their ejections. 

"Considering the uncertainty of to. 

''There is not any optical observation in the B-band available at this epoch. 



Table 2. Parameters of the precession model for the parsec-jet of 3C 120. 



p (yr) a ^ 


an 




n 


10.1 ± 0.8*' 12.5 ± 0.5 


1.3 ± 0.5 


2.6 ± 0.5 


-96 ± 5 



"Measured in the framework fixed at the observer. 

''The quoted error refers to that given by Zhang, Xie & Bai 
(1998). 
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Table 3. Optical flares and radio evolution of the superluminal components of 3C 345. 



Component 


Sop (mjy) 


ATmax (yr) 


<5o 




6^ 


AT43^ (yr)- 


S'op (MJy) 


C3 


1.69 ± 0.12 


3.44 ± 0.27'' 


14.79 


0.50 ± 0.35 


1.22 ± 0.10 


5.6 ± 1.1 


6.0 ± 1.3 


C3a 


2.01 ± 0.32 


1.94 ± 0.37*= 


14.65 


0.76 ± 0.10 


1.31 ± 0.03 


3.7 ± 0.7 


7.4 ± 1.2 


C4 


1.76 ± 0.60 


2.32 ± 0.18= 


18.72 


0.55 ± 0.35 


1.24 ± 0.15 


4.1 ± 1.2 


2.1 ± 0.7 


C5 


2.38 ± 0.22 


1.46 ± 0.16^^ 


22.86 


0.80 ± 0.15 


1.32 ± 0.09 


3.1 ± 0.6 


1.1 ± 0.1 


C7 


2.44 ± 0.38 


1.44 ± O.IQ'^ 


21.35 


0.61 ± 0.16 


1.26 ± 0.06 


2.7 ± 0.5 


1.6 ± 0.2 



Note. — Sopt is the flux density of the B-band optical flares. ATmax is the interval between formation of the 
jet component and its maximum emission at a given radio frequency (see text). So is the Doppler factor at the 
emergence epoch, ar is the spectral index at radio wavelengths and parameter b, defined by equation (16). Finally, 
A^max s-nd ■S'3 are respectively the interval between formation of the jet component and its maximum emission 
and the B-band flux density of the flares in the comoving frame. 

''Considering a = 1; 
t-At 10.7 GHz; 
'^At 22 GHz. 



Table 4. Parameters of a possible black hole binary system in the inner parts of 3C 345. 



Ppt= (yr) 


Pp. (yr) 


r-ps (cm) 


Mp (Mo)'' 


Ms (Mo) b 


4.0 


2.5 


5.5 X 10^'' 


5.0 X 10® 


3.0 X 10® 


5.2 


3.3 


6.6 X 1016 


4.4 X 10® 


3.6 X 10® 


6.1 


3.8 


7.3 X 10" 


4.0 X 10® 


4.0 X 10® 


8.0 


5.0 


8.8 X 10^6 


3.2 X 10® 


4.8 X 10® 


10.0 


6.3 


1.0 X 10^'' 


2.6 X 10® 


5.4 X 10® 


15.0 


9.4 


1.3 X 10^'^ 


1.6 X 10® 


6.4 X 10® 


20.0 


12.6 


1.6 X IQi^ 


1.0 X 10® 


7.0 X 10® 


30.0 


18.8 


2.1 X 10^'^ 


5.6 X 10* 


7.4 X 10® 



-Upper limit; 
''Lower limit. 



